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Abstract 

We investigate the pair production of S-wave heavy quarkonium at the LHC in the color-singlet 
mechanism (CSM) and estimate the contribution from the gluon fragmentation process in the 
color-octet mechanism (COM) for comparison. With the matrix elements extracted previously 
in the leading order calculations, the numerical results show that the production rates are quite 
large for the pair production processes at the LHC. The pt distribution of double J ftp production 
in the CSM is dominant over that in the COM when pt is smaller than about 8GeV. For the 
production of double T, the contribution of the COM is always larger than that in the CSM. The 
large differences in the theoretical predictions between the CSM and COM for the pt distributions 
in the large pt region are useful in clarifying the effects of COM on the quarkonium production. We 
also investigate the pair production of S-wave Be and B* mesons, and the measurement of these 
processes is useful to test the CSM and extract the LDMEs for the Be and B* mesons. 

PACS numbers: 12.38.Bx, 13.85.Ni, 13.60.Le, 14.40.Gx, 14.40.Nd, 14.40.Lb 
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I. INTRODUCTION 



Heavy quarkonium provides a ideal system to investigate both the perturbative and 
non-perturbative aspects of quantum chromodjTiamics. Conventionally, the color-singlet 
mechanism (CSM) is used to describe the decay and production of heavy quarkonium [ij. 
In the CSM, the processes are factorized into two steps. Firstly the heavy quark pair are 
created perturbatively at short distances with the same color and angular momentum as 
the final quarkonium state, and then evolve into the quarkonium non-perturbatively at 
long-distances. There have been, however, some problems in the CSM, e.g., the infrared 
divergences in the calculation of the decay of P-wave quanrkonium and the higher-order 
correction calculation of S-wave quanrkonium [sl, and the surplus J /if) production j4| of which 
the rate is much higher than that of the color-singlet prediction at the Tevatron. The non- 
relativistic quantum chromodynamics (NRQCD) factorization formalism^], which was put 
forward by Bodwin, Braaten, and Lepage, overcame the infrared divergence difficulties in 
the color-singlet modelQ], and gave the proper prediction for the charmonium production at 
the Tevatron j?!. In NRQCD, the heavy quark pair at short distances are not necessarily in 
the color-singlet state but can be in the states with different color and angular-momentum 
from that of the final state quarkonium. The color-octet pair can evolve to the color-singlet 
charmonium by emitting soft gluons. This is called the color-octet mechanism (COM). 

Lots of work have been done to investigate the validity and limitation of the NRQCD 
formulism in heavy quarkonium production. The current experimental results on J / ip photo- 
production at HERA are fairly well described by the NLO color singlet piece except the J/ ijj 



polarizations 



fly, 



lOj. The DELPHI data favor the NRQCD color-octet mechanism for J/ijj 



production 77 —>■ J/ ipX [ill , [l^ . The observed large cross sections of inclusive charmonium 
production at the Tevatron once gave strong support to the color-octet gluon fragmentation 
in NRQCD, but recently it is found that the NLO results in the color-singlet piece can b ring 
an order of magnitude enhancement to the J/ip production rate in the large pt region 13 1 
with longitudinally polarized J/^[l^. The theoretical prediction for pt distribution of the T 
production can properly describe the Tevatron data by including the contributions from the 
NLO results and the real correction part at the next-to-next-to-leading-order in the CSM 
isl . The cross sections of J ftp exclusive and inclusive production in e~^e~ annihilation at B 
factories 16| are much larger than the LO NRQCD predictions 17|, llSl], but the discrepancies 
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seem to be resolved by considering the higher order effects: NLO QCD correctio ns |l9l . |20|. |21| 
and relativistic corrections [22] without invoking the color-octet contributions |20| (discus- 
sions in the light-cone approach can be seen in [23|). Recent developments and related topics 
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in quarkonium production can be found in Refs. 

The above mentioned developments in heavy quarkonium production indicate that the 
situation is far from being conclusive, and further tests for the color-singlet and color- 
octet mechanisms in NRQCD are still needed to clarify various problems involved in heavy 
quarkonium production. 

In order to investigate the effects of the COM on the production of heavy quarkonium, 
it is useful to study processes which heavily depends on the production mechanism. The 
pair production of heavy quarkonium can serve as the desired process. In NRQCD, the 
gluon fragmentation gives the main contribution to the pair production of quarkonium in 
the large pt region of the heavy quarkonium. In the pair production processes, there appear 
two long-distance matrix elements (LDMEs). So the difference of theoretical predictions 
between the CSM and COM could be more obvious. Moreover, because of charge-parity C 
conservation the gluon fusion processes g + g J/ip + Xc and g + g J/ip + r]c are forbidden 
in the CSM. But the gluon fragmentation in the COM can produce these associated final 
states. So to detect two final heavy quarkonium states with different C-parity may give a 
good way to test the COM. 

The pair production of heavy quarkonium at hadron colliders has been studied by many 
authors. The color-octet gluon fragmentation into double charmonium at the Tevatron in 
NRQCD were considered as evidence for the COM . The CSM prediction on the double 
charmonium production was made and it was found that the contribution with pt < AGeV 
in the CSM is dominant [28'] . Only in the large pt region, the CSM and COM give manifestly 
different predictions for the double charmonium production. The Large Hadron Collider 
(LHC) is expected to produce a huge number of heavy quarkonium. Therefore, it is natural 
to investigate the pair production of heavy quarkonium at the LHC. 

At the LHC, it is also interesting to study the production of the double heavy flavored 
mesons Br and B*. The Bc{*) meson production in hadron collisions has been studied in 



QCD 



29 



collisions 
the LHC. 



3C, |3l|]- Some study of i?c(*) pah production was also performed in pp and 77 
32j. It is useful to extend the study to the pair production of the -Bc(*) mesons at 
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In this paper, we study the pair productions of heavy quarkonium at the Tevatron and 
LHC, including J rjc, T and rjb in the CSM. For comparison, the color-octet contributions 
to the pair production of J/ip and T are estimated by considering the gluon fragmentation 
process. We also investigate the pair productions of S-wave and B* mesons where there 
is no contribution from the gluon fragmentation process and the COM contributions are 
suppressed by the small (the relative velocity between quark and anti-quark). Therefore, 
these processes can give a better test of CSM and be used to extract the LDMEs of the Be 
and B* mesons. 

The outline of our paper is as follows. In section 2, some of the definitions and formulas 
are given for deriving the cross sections of the processes. Then the numerical results are 
presented in section 3. Finally, in section 4 we give the summary. 

II. THE FORMULATIONS 
A. Color-singlet part 

At hadron collider, the pair production of heavy quarkonium at the leading-order (LO) 
in the CSM have two subprocesses g + g Qi + Q2 and q + q ^ Qi + Q2- But we just 
consider the gluon fusion process in the calculation since it is the dominant one. There are 
31 Feynman diagrams for the processes g + g —>■ J/tp + J/ip and g + g ^ Bc{B*) + Bc{B*). 
The typical Feynman diagrams are presented in Fig. [H For the process g + g f]c + Vc, there 
are additional 8 Feynman diagrams which are showed in Fig. [2l The Feynman diagrams for 
the process of T and rif, production are as same as the corresponding process of charmonium 
production. 

Following the color-singlet factorization formalism, the amplitude of the pair production 
of S-wave heavy quarkonium is written as 



where the is the spin of the heavy quark in the meson; Ri{0) is the wave function at origin 



M{a + b^ Q1 + Q2) 




xM{a + 6 QiQjiPi = 0; Si, Sa) + QkQi{P2 = 0; S3, S4)) 



(1) 



of the heavy quarkonium; 




is the color project operator; A^j(A|si, S2) is the spin project 
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FIG. 1: The typical Feynman diagrams for g + g ^ J/ip + J/ip- The others can be obtained by 
reversing the fermion Unes or interchanging the initial gluons. As for g + g ^ Bc{B*) + Bc{B*), 
there are same diagrams. 




FIG. 2: The additional typical Feynman diagrams of g + g ^ ijc + rjc- The others can be obtained 
by reversing the fermion lines or interchanging the initial gluons. 



operator as following 



N{X\s,,S2) 



(2) 



where Mqq is the mass of the heavy quarkonium. 

We analytically calculate the amplitude square of these subprocesses and present the 
analytical formulas of the pair production of ric{rib) and J/ip{T) in appendix. For B^. and 
5*, the analytical formulas are too tedious to be presented in this paper. The formula for 
da/dt on the double J/ip production is the same as that in reference 281] . The final result 
can be obtained by convoluting the parton level cross section with the parton distribution 
function fg/p{x) as following 

da{p + p{p) ^ Q1Q2 + X) = j dxidx2fg^/pixi, Hf)fg^/p(ji){x2, fXf)da{gi +92^ Q1Q2, /Ur)(3) 
where the /x^ and /xj are the renormalization and factorization scale. 
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B. Color-octet Part 



As a comparison, we naively estimate the pair production of J /ip and T in the COM by 



using the similar way as in Ref. 27| in which the evolution of the fragmentation function was 
ignored. In the COM, a gluon can fragment into a cc pair with the quantum number ^S^^ 
and then hadronize into J/ip. This process will give large contribution to the double heavy 
quarkonium production in the large pt region and is expressed as 

d^Qi+Q2= dzi dz2Dg^Q^{zi,mQ^)Dg^Q,^{z2,mQ^)dagg{Ei/zi,E2/z2), (4) 
Jo Jo 

where agg is the cross section of two real gluon production; D is the fragmentation function 
for a gluon to fragment into a quarkonium. In NRQCD, this fragmentation function is 
written asjsl 

Dg^Qiz,^,') = J2d,-,n{z,f^'){0^). (5) 

n 

The short distance coefficient can be calculated perturbatively and the result of the LO 
calculation is 

7ras{2mQ) 

The contribution from gg{qq) — * gg subprocesses is calculated and the contributions from 
the feeddown of ip', Xcj(l-P), T(2S') and Xfej(l-P) are also included. The final result is 
expressed as 

rf^Q.+Q. = dagg{^^^^^n{Oi/^eS,)) + {Ot'CS,))Br{4^' - J/^) 

2 

j=o 

Here it is noteworthy that the identical particle factor "2" has been put in the calculation 
of d&gg. 



III. THE NUMERICAL RESULTS AND CONCLUSIONS 

In calculating the numerical results, we choose the following parameters Mc = 1.5GeV and 



Mh = 4.9GeV, and set the renormalization and factorization scale as = /i/ = \r^^Q + Pt- 



For the gluon fragmentation process, the renormahzation and factorization scale are chosen 
as the transverse momentum pt{g) of the gluon with Pt{J/ip) — Pt{g)- The parton dis- 
tribution of CTEQ6L1 331] is used. Therefore, the running of is evaluated by the LO 
formula of CTEQ6. The center-of-mass energies of the Tevatron and LHC are 1.96TeV and 
14TeV respectively. The pseudorapidity cuts on the final quarkonium states are chosen as 
—0.6 < 7] < 0.6 at the Tevatron and —2.4 < t] < 2.4 at the LHC. We use the wave functions 



at origin that are calculated by using the logarithmic potential 



34 



351 ] with quark masses 



almost the same as what we use. The values for them are listed as: 

l^(0)lL-(l5)=0.815GeV^ 
]/^(0)]L(l5) = l•508GeV^ 

l^(0)l'-,(i5)=4.916GeV3. (8) 
In order to give the results of gluon fragmentation processes for comparison, the following 



color-octet matrix elements are used as the input parameters 36|, 1371] 

{Oi^'^CSi)) = 0.39 X 10-2GeV^ {Oj^''\''Si)) = 15 x IQ-^GeV^ 
{OfCSi)) = 0.37 X 10-'GeV^ {Oj^^^\^Si)) = 4.5 x 10-'GeV^ 
^QXco(iP)^3^^)^ = 0.19 X 10-'GeV^ (Of ^'''^(^Si)) = 4.0 x lO-'GeVl (9) 

These LDMEs are extracted from the matching between the LO NRQCD predictions and 
the Tevatron data by using the CTEQ5L parton distribution function with LO a.^ running. 
From the LO formula of ag running of CTEQ6, the corresponding as in the fragmentation 
function for the J/ip and T are chosen as as(Mj/^)=0.286, Q!s(Mx)=0.201 respectively. The 



branching ratios in Eq. ([7]) are taken from the PDG08|38]. 

In the Table m we give the cross sections of pair production of J/ip (T) and Be (B*) at 
the Tevatron and LHC with pt > 3GeV in the CSM. From the table, the cross section of 
each process is enhanced by an order or more in magnitude at the LHC than that at the 
Tevatron. Therefore, the LHC will be a good place to study the pair production processes 
carefully. 

Fig. [3] shows the pt distribution of the pair production for J/ ip and rj^ at the Tevatron 
and LHC. The result of the gluon fragmentation process in the COM is also plotted in the 
figure. We can see that whether at the Tevatron or at the LHC the pt distributions of J/ijj 
and rjc are similar and the numerical results at the LHC are enhanced by an order or more 



TABLE I: The cross sections of pair production of J/V'j Be and B* at the Tevatron and LHC 
with pt >3GeV. 



Final States 


<^Tevatron[nb] 


cTLHcinb] 


VcVc 


3.32 X 10-^ 


2.73 




5.63 X 10-2 


2.83 


VbVb 


1.87 X 10"^ 


7.36 X 10-3 




1.23 X 10-^ 


1.51 X 10-2 


BcBc 


3.86 X 10-=^ 


2.72 X 10-1 


BcB* 


1.00 X 10"^ 


8.37 X 10-2 


B*B* 


8.23 X 10-^ 


7.08 X 10-1 




5 10 15 20 25 30 35 40 45 50 5 1 1 5 20 25 30 35 40 45 50 

P, (GeV) P, (GeV) 



FIG. 3: The pair production of J/ip (solid hue) and rjc (dashed line) in the CSM at the hadron 
colhders. The dotted hne corresponds to the pair production of J/^p that come from the gluon 
fragmentation process in the COM. 
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Tevatron 




5 10 15 20 25 30 35 40 45 50 5 10 15 20 25 30 35 40 45 50 

P, (GcV) P, (GcV) 



FIG. 4: The pair production of T (solid line) and r/b (dashed line) in the CSM at the hadron 
colliders. The dotted line corresponds to the pair production of T that come from the gluon 
fragmentation process in the COM. 

in magnitude at large pt region than that at the Tevatron. Therefor, the LHC will provide 
a chance to measure the pt distribution of J/ip pair production. Comparing the J/ip pair 
production in the CSM with that in the COM, the formal is dominant as the pt is smaller 
than about 8GeV. And the result in the COM become dominant and even larger than that 
in the CSM for three orders in the large pt region. 

The Pt distributions of T and rjf, production are shown in Fig. HI The pt distribution is 
enhanced by an order or more in magnitude at the LHC than that at the Tevatron. But 
unlike the case of J/ ip pair production, the pair production of T in the COM dominate over 
that in the CSM in the whole pt region, and even is two order or more in magnitude larger 
than that in the CSM at large pt. 

Because the Be and B* are consist of quarks with different flavors, there is no contribution 
from gluon fragmentation processes in the COM. So in Fig. El we only give the pt distribution 
of the pair production for B^Bc, B^B* and B*B* in the CSM. The pair production of B*B* 
is dominant in the whole pt region. There are little difference between the production of the 
three final states with double-heavy flavor mesons in the CSM. 

From the Table I and the figures, it can be seen that all of the cross sections are enhanced 
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5 10 15 20 25 30 35 40 45 50 5 1 1 5 20 25 30 35 40 45 50 

P, (GeV) P, (GeV) 



FIG. 5: The production of B^Bc (solid line), B^B* (dashed line) and B*B* (dotted line) in the 
GSM at the Tevatron and LHG. 

when the center-of-mass energy is increased. This is because that with the fixed pt, the larger 
the is, the smaller the momentum fraction x of the parton is. In small x region, the 
parton distribution function of gluon increases rapidly. 

IV. SUMMARY 

In this paper, we have investigated the leading order pair production of S-wave heavy 
quarkonium at hadron colliders in the color-singlet mechanism (CSM) and estimated the 
contributions from the gluon fragmentation process in the color-octet mechanism (COM) 
for comparison. With the matrix elements extracted previously in leading order calculations, 
the numerical results show that the production rates are quite large for the pair production 
processes at the LHC. The pt distribution of double J/t/j production in the CSM is dominant 
over that in the COM when pt is smaller than about 8 GeV. For the production of double 
T, the contribution of the COM is always larger than that in the CSM. There are large 
differences in the theoretical predictions between the CSM and COM for the pt distributions 
in the large pt region, and this is useful in clarifying the effects of COM on the quarkonium 
production. Furthermore, since to produce a pair of quarkoniums with different C-parity 
is forbidden in the CSM at the leading-order, the observation of these processes could be a 
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positive support for the COM. We also investigate the pair productions of S-wave Be and 
B* mesons, and the measurement of these processes may be useful to test the CSM and 
extract the LDMEs for the B^ and B* mesons. 

39|, a paper appeared in which Qiao, Sun, and Sun 



After our work was completed 



calculated the double J/ip production at the LHC They focused on the polarizations of the 
double J/ip. We focused on the cross sections of double heavy quarkonia J/ip, Vc, T", Vb, as 
well as the double heavy flavored B^, and B* mesons. Both the two papers discuss the test 
of the COM. Our color octet and color singlet double J/ip cross sections are consistent with 



their result [39 1. 
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Appendix 

In this appendix we give the parton level cross section for the pair production of ric{rib) 
and J/iIj{T) respectively. Here m is the mass of the corresponding meson and the s, t, u are 
the Mandelstam variables defined as 

s = {ki+ /c2)^ 
t = ih- p)\ 

u = ih- P)\ (10) 

where fci, k2 and P are the momenta of the initial gluons and one of the final quarkonium 
states. 
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^- 9 + 9 ^ Vc + Vc 



162m2s8(m2 - t)H\-2'm? + s + ty{-'m? + s + tY 
+m)m^^ + 256i(27s^ + 14257^5^ + 81792^^5 + 101412i^)m2^ - 64i(543s^ + 72538is^ 

+687780^^5^ + 1876968t^s + 1498176t^)m22 + 4(9s^ + 18144ts''' + 18144ts^^ 
+1018456^2^4 + I3409824i3s3 + 60363088^^5^ + 104825088i^s + 60673536t^)m2° 
-4(45s^ + 20364^5^ + 693676^^55 ^ I0769920t^s^ + 70075680^^^ + 199182848t^s2 
+248330880tls + 111310848t^)mi^ + (369s^ + 53568t/ + 1532448t2/ + 24574336f'^s^ 
+210621248tV + 870605312^^5^ + 1761630976^^5^ + 1686196224i^s + 610384896t^)m^^ 
-4(995^^ + 5496^5^ + 172056^^5^ + 2626424i^s^ + 27129120iS^ + 152150816i^s^ 
+446139360t^s^ + 687430400t^s2 ^ 527406336t^s + 158754816t^)mi^ + 2(117s^° 
+3420fe^ + 131584^^58 + 1843176i^s^ + 20064012i^s^ + 143086496i^s^ + 571404544i^s^ 
+1267269888^^5^ + 1553349056^^5^ ^ 984517632^*^5 + 251817984i^°)m^2 _ 4(18511 
+588ts^° + 21495t2s9 + 304486t^s^ + 2927138t^s^ + 23491008t^s® + 120913176t^s^^ 
+363365760iV + 636265120^^5^ + 639890816i^52 + 342648576i^°5 + 75727872i^^)m^° 
+(95^2 ^ 810^5^^ + 20676i2sio + 360708^^5*^ + 3111538^^^ + 23098608^^5^ 
+138861808t^s^ + 537887488t^s^ + 1266810688t^s^ + 1806816000t^s^ + 1525860352t^°52 
+702443520^^^5 + 135945216i^2)m^ - 4^(455^^ ^ 903^^11 + 18636^^510 + 181253^^5*^ 
+1173996iS^ + 7001216^^5^ + 32407688^^5^ + 98050752i^5^ + 186725360i^5^ 
+222162048tV + 160423104^^^ + 64384128tiis + 11031552ti2)^6 ^ 2t{'ds^^ + 396^^2 
+4962i2sii + 58165^^51° + 384536i^5^ + 2056988^^5^ + 10083672iS^ + 37041746^^5^ 
+90302176^^5^ + 142532816^*^5^ + 144160704fi°5=^ + 90317984^^^52 + 31965696i^2g 
+4893696t^^)m^ - 4t2(45si^ + 'ims^^ + 27m^s^^ + 20392f'^si° + 103737t'^s^ 
+480982^^5^ + 1974841t^5^ + 5914228^^5^ + 11914204^^5^ + 15898752^^5^ 
+13880896ii°5=^ + 7636608^^^52 ^ 2406528^^2^ ^ ^'^nm'^^)rn^ + ^2(5 + tf{ld>s^'^ 
+90ts" + 447t2sio + 3162t^s^ + 14485t^s^ + 57520t^s^ + 241296t'^6'^ + 755200tV 
+1481344^^5^ + 1774080^^5^ + 1267200ti°52 + 497664^^^5 + 82944^^2)) (i^) 
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? = . '^^'"'!^^°^r ^(7776m-- 432(73. + 216t)m- 

+6(90855^ + 60336^5 + %hhZU^)m^^ - 16(36295^ + 37686^8^ + 117855^^5 
+106920i^)m^^ + 4(11927s^ + 151588^5^ + 745674^^5^ + 1470960i^s 
+962280t^)m^^ - 4(7761s^ + 109608ts^ + 699467t2s^ + 2173908t^s2 
+3055320^^5 + 1539648t^)m^^ + 2(6952s^ + 117893is^ + 897043i^s^ 
+3741980t^s^ + 8278410^^5^ + 8872416i^s + 3592512i6)m^2 - 2(1899s^ 
+43398ts^ + 405618t^s^ + 2113568t^s^ + 6394090t^s-^ + 10762584^^5^ 
+9189936^^5 + 3079296i^)m^° + (587s® + 19710is^ + 244772^^5^ 
+1603468t^s^ + 6229962iV + 14478304i^s^ + 19359816i^s^ + 13582080i^s 
+3849120t®)m® - 2(20s^ + 1185ts® + 22153t^s^ + 193780t^s^ + 965358tV 
+2928368i^s^ + 5431786^^5^ ^ 5949528^^5^ ^ 3508920i®s + 855360i^)m^ 
+(s^° + ms^ + 3756^^5® + 52062^^5^ + 353472tS^ + 1398834i^s^ 
+3421754t'^s^ + 5210968t^s^ + 4784622t^s2 + 2414880t^s + 513216t^°)m^ 
-M^{s + i)2(9s^ + 649ts^ + ^A&^i^s^ + 29630i^/ + 74435i^s2 + 105156^^5^ 
+77868iS + 23328i^)m^ + 2t\s + i)^(349s^ + 2304is^ + 6192^^52 
+7776t^s + 3888t'^)) (12) 
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